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'. Abstract -In pseudo integrable systems diffractive scattering caused by wedges and im- 



On 
^— > 



purities can be described within the framework of Geometric Theory of Diffraction (GDT) 
in a way similar to the one used in the Periodic Orbit Theory of Diffraction (POTD). We 



'"O ! derive formulas expressing the reflection and transition matrix elements for one and many 

s ■ 

O ' diffractive points and apply it for impurity and wedge diffraction. Diffraction can cause 



backscattering in situations, where usual semiclassical backscattering is absent causing 



j> ', an erodation of ideal conductance steps. The length of diffractive periodic orbits and 
00 ■ 

00 diffractive loops can be detected in the power spectrum of the reflection matrix elements. 

■ The tail of the power spectrum shows ~ l// 1 / 2 decay due to impurity scattering and 

O ; 

O" < ~ 1/ 1 3 / 2 decay due to wedge scattering. We think this is a universal sign of the presence 



of diffractive scattering in pseudo integrable waveguides. 
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£H ■ In recent years, semiclassical methods became very popular in describing devices operating in the 

O '. 

Qmesoscopic regime. There are two entirely different sets of theoretical tools which have been used in 
. £naj wide range of applications. One of them is a WKB based short wavelength description [[l], ^ where 

^classical trajectories, chaos, regularity and analytic properties of the potential play a major role. The 

d | 

other is based on random matrix models or on averaging over random Gaussian potentials . The first 
approach is designed to describe clean systems where the potential depends smoothly on coordinates and 
parameters, while the second assumes a system densely packed with impurities causing wild fluctuations 
of the potential on all length scales. 

Fortunately, advances in manufacturing and material design reduced the average number of impurities 
exponentially since the eighties and this trend is expected to continue in the future. Accordingly, a 
realistic semiclassical theory should be capable to treat systems with low number of impurities and not 
only completely clean or completely dirty ones. Another demand is that human designed structures, 



ay 




Figure 1: Electron reservoirs (Ohmic contacts) connected trough strip waveguides, a, Ideal strip waveguide 
with parallel walls, b, Ideal strip waveguide with parallel walls and one pointlike impurity, c, Broken-strip 
waveguide. 

straight sides and eventual wedges should be a natural part of the description. These systems, neither 
chaotic nor regular, form the intermediate class of pseudo integrable systems. In the present paper we 
will demonstrate that diffractive scattering caused by wedges and impurities can be described within the 
framework of Geometric Theory of Diffraction [|], [5], || [7|, H (GDT) in a way similar to the one used in 
the Periodic Orbit Theory of Diffraction || [H], [ll| (POTD). The POTD is an extension of Gutzwiller's 
periodic orbit theory for non-smooth systems with diffractive scattering which has been used successfully 
during the last three years. 



1 Examples of pseudo integrable mesoscopic devices 

The simplest mesoscopic device is a two dimensional wave guide strip of electron gas, formed on a GaAs 
heterostructure |T^], connecting two electron reservoirs (Fig. la). This system can be considered also 



as a paradigm of a more sophisticated two or three dimensional clean quantum wire ^5) made of other 
materials, for example carbon based structures with delocalized electrons. Here we will consider two 
slightly modified versions of this simple quantum wire, the quantum wire with one point like impurity 
(see Fig. lb) |16|], the broken-strip waveguide (see Fig. lc) fT7|| . Our motivation in selecting these 



systems for investigations is that any long wire with straight sides, many wedges and impurities can be 
built up as a sequence of these building blocks. For simplicity, the walls of the wave guides are assumed 
to be completely hard and the wave function vanishes on them. The potential is assumed to be zero 
inside the waveguide. 

2 The Landauer Formula 

The conductance of the non-degenerate electron gas in a waveguide considered here is given by the 



Landauer formula [18, 19 1. According to this theory, incoming and outgoing wave functions in the leads 



far from impurities or wedges can be decomposed into incoming and outgoing quantum modes 

i() n {x,y) = <5> n (y)e lknX , 

®n{y) = y^sin(7rm//W0, (1) 



where k n = ^2mEp — (fiim/W) 2 /% is the wave number of the propagating planewave, W is the width of 
the lead and E F is the Fermi energy. If the Fermi energy is less than E n = ^(h%n/W) 2 the wavenumber 
k n becomes imaginary and the n th channel becomes closed preventing wave propagation. 
The Landauer formula for the conductance G(Ep) at Fermi energy Ep is then 

2e 2 2 

G(E F ) = —— ^ \tmn\ j (2) 
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where t nm is the transition probability amplitude from the incoming cannel n on the entrance side to the 
outgoing channel m on the exit side. In case of infinitely long leads summation goes for open channels 
of both sides only. The transition probability is given by the projection of the Green function over the 
transverse wave functions & n {y) on the entrance lead for the incoming modes and <&' m (y') on the exit 
lead for outgoing modes 

t nm = 2i{k n k' m ) 1 ' 2 j dydy'$ n (y)G + (x,y\x',y')& m (y'), (3) 



and x and x' lie anywhere on the entrance side and exit side, respectively (See Fig. la-c). (From here 
on we will use units such H = 2m = 1.) Reflection between two modes n and m on the same side is 
given by 

r nm = 5 nm - 2i(k n k m ) 1/2 J dydy'® n (y)G + (x, y\x, y')<$> m (y'), (4) 

where x lies anywhere on the entrance side. For open channels n of the entrance side the transmission 
and reflection amplitudes fulfill the sum rules 

\tnm\ ~\~ l^nm'l = lj (5) 



where the two summation goes for channels m on the exit side and w! on the entrance side. These are 
the consequences of the probability and current conservation. Thanks to these relations, in infinitely 
long wave guides, the conductance can be expressed with the reflection coefficients 

G(E F ) = ^(N-Y / Km\ 2 ), (6) 

where N is the number of open channels on the exit side. 

The direct connection between the Green function and the conductance makes it possible to develop 
a semiclassical approximation for the conductance based on the semiclassical Van Vleck - Gutzwiller type 
expressions for the Green function. Next, we summarize the main tools of the semiclassical description. 

3 The semiclassical Green function 

The Schrodinger equation in the waveguides leads to the Helmholtz equation 

(A + k 2 )^(q) = 0, (7) 



where the notation k = has been introduced. The wavefunction vanishes on the boundary as 

we discussed above. The energy domain Green function of the Helmholtz equation is defined by 

(A + k 2 )G(q,q>,E) = 5(q-q>), (8) 

and G(q, q', E) should vanish on the boundary of the waveguide too. In absence of walls in two dimensions 
the outgoing Green function is given by 

G+(q,q>,E) = - t -Hl ) 1 \kd(q,q')), (9) 

where Hq 1 \x) is the Hankel function of first kind and d(q, q') is the distance between q and q' . The semi- 
classical (H — > 0) approximation of this free space Green function can be recovered from the asymptotic 
form of the Hankel function for large argument 

G + (q, q', E) = 1 e ikd( q , q ')-i3n/4_ (1Q) 

'87rkd(q,q') 



A very useful optical interpretation of this formula can be given by tracing the ray connecting q with q' . 
The phase of the Green function is the classical action calculated along the ray 



f p(q"W = kd(q,q>) 

Jq 



The amplitude of the Green function is the square root of the intensity 



I(r) = —L- (11) 
V ; 8nkr K ' 



of a radiating point source where the distance from the source is r = d(q,q'). The phase factor e~* 37r / 4 
is the Maslov index of the caustics singularity standing right in the source point q. 

In the presence of walls, in general, more than one ray from q can reach the point in q'. Then the 
semiclassical Green function is a sum for the raywise Green functions calculated along the rays 

G(q,q',E)= £ G r (q,q',E). (12) 

Vrq— *q' 

The phase of the Green functions G r is given by the classical action plus the Maslov index including the 
caustics in the source point. The amplitude is the square root of the intensity observable in q' coming 
form q along the ray. In the endpoint we will see the mirror image of the source trough perfect mirrors 
given by the walls. The intensity in this case is given by the same formula (0), but now the effective 
radius of the spherical wave coming from q should be used. If the walls are non curved this radius is 
simply the distance d r (q, q') of q and q' along the ray. In this case the Green function is given by 

G r Q (q, q', E) = 1 ' e ikd r (q,q')- l3 n/^ ^ 

'8nkd r (q, q') 



where n r is the number of bounces along the ray. For the semiclassical expression of the Green function 
in presence of curved walls the effective radius is given (Tl| by the product of the distance between the 
source and the first bounce point l\ times the stretching factor A 



eff 



hA. (14) 



The stretching factor is the analog of the optical magnification factor of curved mirrors. It can be 
calculated as 

A = n f[(l + 4k), (15) 

i=l 

where Zj is the length of the free flight between the i th and the (i + l) th bounces and k is the Bunimovich- 
Sinai curvature. The Bunimovich- Sinai curvature is defined recursively 

4 = kT + 2/rjCos(0j), 

= 4/(1 + 4k), (16) 

where Ti is the radius of curvature of the wall in the point of bounce and 0, is the angle of incidence. The 
initial condition of the recursion is rC[ — 1/h, which is the Bunimovich-Sinai curvature of the wavefront 
arriving at the first bounce point from the source. The Green function reads as 

G5(?, q', E) = ^(^W4_ (17) 

Vo7TK<iA r 

In the special case of non-curved walls the stretching factor (optical magnification) is A = d r (q,q')/li 
yielding the expression 



4 Geometric Theory of Diffraction 



The semiclassical approximation works well when changes of the potential are on a length scale much 
larger than the wavelength of the electron 1/k |[20|| . Even if the manufactured system in large satisfies this 
condition, there are isolated points such as corners or impurities where this cannot be the case. These 
points require special treatment || F2], [TB|]. They break the smooth wavefronts of the semiclassical 



wave propagation and create diffracted waves coiling out from them. Such diffractive points can be 
considered as new wave sources whose strength is proportional with the strength of the incident wave. 
Their contribution of such a specific ray to the Green function is given by 

G d (q, q', k) = G (q, q , k)V(4> in , 4> out , k)G (q , q', k), (18) 

where G (q, q , k) is the Green function calculated along the ray connecting the source point x and the 
diffractive point qo, T>((j)^ n ,(j) r out ) is the diffraction constant which can depend on the incoming and the 
outgoing angle of the ray and the energy. Go((Zqj q '■> k) is the Green function calculated from the source 
to the point of observation q'. If there is more than one ray connecting q with go an d qo with q' then 
each ray configuration from q to q' will contribute to the Green function according to fll^D- 

The diffraction constant can be determined for different geometric situations. For wedge diffraction 

1 1 



v( j. x \ sin(vr/n) 



'ini 'rout j 
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cos(7r/n) - cos((0 in - (fi out )/n) cos(7r/n) - cos((7r + 4> in + <p out )/n) 



(19) 



where (2 — n)ir is the opening angle of the wedge. For details we refer to Ref. |5|, |TT], For diffraction 
on impurities the s wave scattering is dominant and the diffraction constant D is isotropic, depending 
only on the energy. Typical example is a ^-shaped potential with scattering strength T>. 

5 Semiclassical transmission and reflection 

Semiclassically the transmission (|3|) and reflection (||) matrix elements are calculated by replacing the 
Green function in (0) and (f|) with the semiclassical expressions (|17|) or flT3|) . The wave functions $ n (y) 
and & m (y') can be decomposed to exponentials 



*M) = I^/A( e -V/^_ e -W/^. (20) 
The transmission and reflection matrix elements are then sums of four subintegrals 

nm nm ' nm ' nm ' nm ' 

trim — t nm + t nm + t nm + t nm , (21) 



where 



nm 



5nm - iOW™) 1 ^ / dydf/e^l w ^l w -^l 2 G + {x, y\x, y'), (22) 
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where z/ = ±1 ± 1. By using (]i~7f) or (p~3|) in semi classical approximation one can calculate these integrals 
with saddle point approximation. The saddle point conditions for various t^-s are 

kd y d r (x,y\x' ,y') = ^fnn/W and 

kd y id r (x,y\x' ,y') = ^fnm/W', (23) 

and for r^-s are 

kd y d r (x,y\x,y') = ^fnn/W and 

kd y 'd r (x,y\x,y') = ^pirm/W. (24) 

If we introduce the initial 0j„ and final </>/j n angle of the rays with respect to the direction of the leads 
these conditions are equivalent with selecting rays whose initial angle is determined by 

sin(0 m ) = (25) 

and their final angle is 



sm(0 /in ) = =F ttjt; for *■ 



nmi 



kW 

sin((j) fin ) = =Ftt77 for r nm, (26) 
kW 

respectively. 

For irregularly shaped wave guides one can in general find rays starting and ending with given initial 
and final angles. We can call such systems "chaotic". In wave guides formed by regularly shaped walls 
it is possible that in the bouncing process a components of the momentum is conserved in addition to 
the energy conservation. In these systems only one of the saddle point conditions can be satisfied. We 
can call them "integrable". One example is the waveguide on Fig. la. There are systems, where chaotic 
and integrable behaviour coexist. For certain initial angle one can find rays to any final angle, while for 
other initial conditions not. These systems are "mixed". The conductance of chaotic, integrable and 
mixed systems are relatively well understood. 

The simple, typically human designed systems of Fig. lb and Fig. lc do not belong to any of these 
categories. We may call them "pseudo integrable" systems. These systems are not chaotic, since there 
exist is no ray in general connecting a given initial and a final angle. An initial angle allows only a finite 
number of possible final angles. So, they are more like integrable systems in this respect. However, when 



we take into account the diffractive points, rays can leave them in any direction. So, we can always find 
diffractive rays connecting any initial angle with any final one. In this respect, these systems are more 
chaotic like. 

In the wave guides on Fig. lb and lc starting from the left with a given initial angle and bouncing 
on the walls, the electron can end up on the right with a definite final angle. Return to the left is 
impossible. Except for a classically "measure zero" fraction of rays hitting the impurity or the wedge. 
The contribution of these rays to the transmission and to the reflection can be computed by substituting 
the diffractive part of the Green function (|18D into (p2|). The length of the ray connecting (x,y) with 
(x',y') or (x,y') via the diffractive point (x ,y ) is 

d(x,y\x',y') = d{x, y\x , y ) + d(x , y \x', y'), 

d(x,y\x,y') = d(x, y\x , y ) + d(x , y \x, y'). (27) 

Since the position of the diffractive point (xo,|/o) is fixed, the saddle point conditions for the initial and 
final segments decouple for t^-s 

kd y d(x,y\xo,y ) = ^irn/W and 
kd y ,d(x ,yo\x',y') = ^irm/W, (28) 



and for the r, 



kd y d(x,y\x ,y ) = ^nn/W and 
kd yl d(x ,yo\x,y') = ^ixmjW. (29) 



i£ = %fe=^ I dydy>e ±m ^ w±m ^/ w '-^ 2 G (x,y\x ,y )V(^ 



The saddle point integrals 

Vww 



r ±± 
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l{knk ™) 1 ' 2 f d^'e^/^™^^ (30) 



then can be carried out separately. The Kronecker delta disappears, as usual |L9], 0, due to the trivial 
direct rays from y to y'. Assuming straight walls we get 

(_-l\n r 

+±± _ _,,'"TV/A rh \ ^ ' ik(d(x,y\x ,y )+d(x ,y \x',y'))±inTTy/W±imTTy'/W'-iuTT/2 

iJWW'knk' 



(_-l\n r 

±± _ j_AT)(sh sh \ : 1 ik(d(x,y\xo,yo)+d(xo,y \x,y'))±in-Ky/W±irmry' /W-ivir/2 



where n r is the total number of bounces on the straight wall and y and y' are the solutions for the saddle 
point conditions ( |2lf ) and (P9"|). The ± signs in y± and y'± have been suppressed for brevity. 



A ray with only one diffractive point is the most elementary possibility for diffraction. More com- 
plicated diffractive rays are those, which return to the diffractive point several times. These events are 
described by the following Green function 

G d (x, y\x', y') = G (x, y\x , y )G (x , y \x', y') (v((f) in , <f) out ) + ^ V{4> in , (f) r )G r ((f) r , 4>' r )V(4>' r , ou t)J , 

(32) 

where 0j n and cj) out denote the angles under which the incoming and outgoing rays reach first and leave 
finally the diffractive point, <fr r and <p' r stand for starting and ending angles of ray loops starting and 
ending in the diffractive points and Go(0 r , 4>' r ) denotes the Green function computed along this ray loop. 
The expression in the bracket can be considered as an effective diffraction constant, renormalized by the 
interaction with the environment. The closed ray loops starting and ending in a diffractive point were 
introduced in Ref. ||] first and are called diffractive periodic orbits. The summation for r is a summation 
for all the possible diffractive periodic orbits starting and ending in the diffractive point 

V(4> in , 4> out ) = V(<f) in , <p out ) + £ A p e ikd », (33) 

v 

where d p is the length of the diffractive periodic orbit and A p = \V((p in , (fi p , k)G r ((f> p , 4>' p )V((p p , 4> ou t)\- The 
transmission and reflection matrix elements then can be written in terms of the renormalized diffraction 
constants 

(_-\\n r 

+±± _ _oT)(^h rh \ 1 l Ak(d(x,y\xo,yo)+d(xo,y \x',y'))±inT7y/W±imT7y' /W-iutt/2 

L nm ~ bL -'\H'%mH'out) / e j 

Ay/WW'k n k' m 

(_-\\n r 

r ±± _ -^ijjffj)^ (j) out ^j c ik(d(x,y\x Q ,y )+d(x ,yo\x,y'))±in-Ky/W±im-Ky' /W-ivn/2 

4Wy/k n k m 

The expression derived here is the main result of the paper and expresses the diffractive part of the 
transmission and reflection matrix elements in terms of diffractive periodic orbits and the rays reaching 
and leaving the impurity. 

In case of more than one diffractive points qo, ...,qk, one can imagine more complicated situations. 
The ray from y can reach qi first, then jumps to q s and goes to y'. Beyond the direct ray between / 
and s, there are rays which start in / and reach s in different ways, possibly trough diffraction on other 
diffractive points or trough bouncing on the walls or trough some combination of these. These type of 
processes are described by the Green function 

G d (x,y\x',y') = ^2G (x,y\x h y l )G (x s ,y s \x',y') ^£>(0^/, (p r )G r ((p r , (p' r )V(4>' r , 4> s ^ y >) \ , (35) 

l,s \ r / 

where G>(0 r , (f>' r ) is the Green function computed along the ray leaving / at angle and reaching s at 
angle (f)' . We can introduce the notations 

Vi^y^i^^y^k) = ^P(0^ / ,0 r )Go(0 r ,0^)P(0^,0 s ^2 / /), (36) 



where the diagonal terms T>u((j)y-*i, <f) s ->y', k) are defined as the renormalized diffraction constants defined 
in (^). Finally, we can compute the transmission and reflection matrix elements as before and get 



+±± _ A T), I A, , A, U\ v / ik(d(x,y\xi,y l )+d(x s ,ys\x',y'))±iniTy/W±innry' /W'-iv%/2 

u 4JwW'k n k' m 
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( 1 \n__ 

±± _ ijSPT),(A , fk \ V > ik(d(x,y\x l ,y l )+d(x s ,y 3 \x\y'))±inwy/W±imi<-y'/W-ii'T!-/2 /o 7 \ 

nm — / , t^ls\9in,h 9out,s) ^ — , c . ^'J 

This is the most general expression, within the framework of the GDT, for the transmission and reflection 
matrix elements. 

Next, we are going to apply the theory developed here for the examples announced in the introduction 
and show which are the consequences of diffractive contributions for the conductance of these systems. 

6 Erosion of conductance steps 

We have calculated the conductance of an infinitely long waveguide with a point like impurity in its 
center and the broken-strip waveguide (see Fig. 2). The details of the calculation will be published 
elsewhere |2"T|. The geometry in both cases is such, that a classical electron starting from the left 



of the device cannot return back via classical bounces on the walls ||23|| . The usual semiclassical theory, 
in the absence of returning rays, predicts no reflection for these systems and the semiclassical reflection 
matrix elements r nm all vanish. According to (^) the conductance of such systems should coincide with 
that of an ideal straight waveguide G(E E ) = ^jr-N. As we can see on Fig. 2, this is not the case. 
The difference between the sharp ideal staircase and the actual conductance is the consequence of the 
non-vanishing reflection. In case of the pointlike impurity the reflection can be determined exactly and 
various terms can be interpreted from a diffractive point of view. The exact reflection matrix elements 
for an impurity in an infinite straight waveguide are 

r nm = -S»(tt))gm(tt>) . L^ e^'"^^! V -, (38) 

2t^k n k m 1 - VG E {x , y \x , y ) 

where the exact Green function of the empty guide 

G E (x, y\x\ y') = £ ^M^l e ^-A (39) 
has been evaluated on the impurity and the singularity has been removed 

G E {x , yo\xo, Vo) = Km (g e {x + e, y \x Q , y ) + ^Y (ke)^j , (40) 



where Yq(x) is the Neumann function. By using (|20"D we can determine the exact elements: 



i D 

^,±± g±inyo/W±imyo/W ^ik \x—xo\—w-k/2 (41) 

4W^/k n k m l-VG E (x ,y \x ,yoY 





Figure 2: Conductance as a function of the Fermi wavelength in units of a, For the infinite strip with 
impurity (T> = 1). b, For the broken-strip waveguide. The angle of the strips measured from the line connecting 
the two wedges were tt/6 and ir/2. 



Figure 3: Typical diffractive rays, a, Primary rays reaching and leaving the diffractive point created by the 
pointlike impurity, b Primary rays reaching and leaving the diffractive point in the wedge, c Primary ray and 
the ray bouncing forth and back between the diffractive points. 

These terms can be interpreted as contributions from various rays reaching the impurity. For example 
a ray starting in x, y under angle sm(<f) in ) = nir/Wk will reach Xq,i/o by flying a distance \x — Xq\ in 
horizontal direction and \uq — y\ + rW in the vertical direction, where r is the number of bounces during 
the flight, the total length of the ray is 

d(x, y\x , y ) = (\x - x \ 2 + (\y - y\ + rW) 2 ) 1/2 , 

and 

cos(0 in ) = \x - x \/d(x,y\x ,yo), sm(4> in ) = (\y - y\ + rW)/d(x, y\x , y ). 
Using the relation k n = k cos(0j n ) the full exponential can be written as 

k n \x — Xq\ + riTcy /W = kd cos 2 '(fan) + k sin(0 in )((isin(0 in ) — rW + y) = kd — mix + niry/W. 

Replacing the phase terms in ( p]) with this and analogous expression for the other ray (see Fig. 3a) we 
recover (^4j) exactly with renormalized diffraction constant 

V 

V = * . (42) 

1 - VG E (x ,y \x ,yo) 

It is very reassuring to see diffraction theory reproducing the exact result in this case. 

The regularized Green function can be written semiclassically as a sum for all trajectories starting 

on and returning to the singularity |24 



G E (q ,q )^ E -V4jr C ■ (43) 



Accordingly the renormalized diffraction constant 

V = V + V 2 G {x , y \x , y ) + V 3 G 2 (x , y \x , y ) + ... (44) 

can be viewed as a sum for all returning rays having one, two, ... etc. diffractions on their way. 

It is very enlighting to study on this simple model how diffraction erodes the conductance steps 
close to Fermi wavelengths where new channels are opened k — > im/W (k n — > 0). Around these Fermi 
wavelengths the regularized Green function is dominated by a single term 

Goto,*)* *-™?"™ . (45) 



2ik 



The squares of the reflection matrix elements can be approximated by 

2 _ k n |$ m (|/o)| 2 1 

The nondiagonal reflection coefficients are small and vanish for k n — > +0, while the diagonal is a 
Lorentzian of width A„ = T>\$ n (y )\ 2 /2 

1 

1 + (k~/A n ) 

Accordingly, the conductance steps have an eroded shape 



I 2 = 1 I n. / a ^ k n > 0- ( 4 7) 



We can see, that the average width of the erosion of the conductance steps is a primarily information 
about the level of diffraction in a waveguide (see Fig. 2a). 

The diffraction theory is not exact for the broken-strip waveguide, but the conductance shows quali- 
tatively the same behaviour (Fig. 2b). The main difference is that the diffraction constants here depends 
on the angle and reflection matrix elements are more complicated. The primary diffractive rays reaching 
and leaving the wedge are depicted on Fig. 3b. For k just above nn/W the initial angle is close to 
(p in »i 7r/2 e.g. the electron is bouncing almost perpendicularly on the walls of the guide and cannot 
reach the upper wedge via bounces. On Fig. 3c a ray reaching the lower wedge by bounces and then 
going to the other wedge by diffraction is shown. When k gets a bit larger, the initial angle decreases 
and rays reaching the other wedge directly should be considered. The sudden appearance of the new ray 
causes an abrupt change in rn at k — 1.5 on Fig. 4. 

For larger k » nn/W Fermi wavelengths the reflection matrix elements decrease rapidly and show 
oscillations small amplitude oscillations dominated by the oscillations coming from individual diffractive 
orbits (see Fig. 4). 





Figure 4: The reflection matrix element r\\ as a function of the Fermi wavelength in units of ^ a, for the 
strip waveguide with impurity b, for the broken-strip waveguide. (Re rn(fc), Im ru(k) are denoted by solid 
and dashed lines, respectively) 



7 Diffractive periods in the power spectrum 

For k » nn/W both k n and k m can be approximated with k. The primary diffractive rays are almost 
parallel with the walls of the guide 0^0. Since the Green functions of the rays have a prefactor 
proportional with l/y/k, multiple diffractions are largely suppressed. The reflection matrix elements 
can be well approximated in leading 1/k order by keeping rays with one and two diffractive points only. 
Rays going trough two impurities while traversing the system are the most trivial examples. The first 
non-trivial contributions come from rays reaching a diffractive point and leaving it after making a loop. 

In case of the point impurity this means that we keep terms proportional with the first and second 
power of the diffraction constant only: 

rum ~ -$n(yo)Myo)^e ifc "l*-* ol+iM*-* ol (v + V 2 G E {x Q ,y Q \x Q ,y )) . (49) 

The regularized Green function can be written as a sum for rays starting and ending on the impurity 
yielding 



Turn ~ -Qniyo^miVo)^ 2 ^-^ [V + V 2 ]T J_L~ ^(^o)-^ (5Q) 

Beyond the trivial phase factor e t2k \ x - x o\ coming from the ray reaching and leaving the impurity, the 
most important oscillations are caused by rays making a single loop starting and ending on the impurity. 
The length of these orbits can be calculated and can be grouped in three subcategories 

d r = 2y + 2Wr, 

d r = 2(W - y ) + 2Wr, (51) 
d r = 2Wr, 

where r = 0, 1, 2, .... The amplitude of the oscillations is proportional with \ j \fd~ r ~ 1/y/r. 

In the broken-strip waveguide the situation is somewhat more complicated. Loops starting and 
ending on the wedges can do it typically by making additional diffraction on the other wedge as shown 
on Fig. 2. But this diffraction process is somewhere in between a usual perpendicular bounce back from 
a wall and a diffraction in the corner. Therefore, the contributions of these orbits is not suppressed after 
one bounce. To see how the amplitude of these orbits decreases, we Fourier transformed the large k tail 
of the reflection matrix element rn and calculated the power spectra 



MO 



oo 

dkrn{k)e- m 

ko 



(52) 



where k w 2ir/W has been chosen. We have multiplied r u with a phase factor e l2kL in order to remove 
the main oscillation coming from the distance L between the lower wedge and the starting point of our 
coordinate system. On Fig. 5 we can see that the contribution of the suspected orbits is the most 
significant. Their amplitude decay as ~ 1/r 3 / 2 , which is faster than it is in the impurity case but slower 
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Figure 5: The power spectrum of the reflection matrix element rn for the broken-strip waveguide. The heights 
of the peaks are proportional with the Z -1 ' 5 . 

than the exponential decay predicted for multiple diffraction. 

We can conclude, that in these two pseudo integrable situations the non-trivial part of the power 
spectra of the reflection matrix elements shows a ~ type scaling with (5 — 1/2 and (3 = 3/2. In a 
similar but chaotic situation the power spectrum decays much faster due to the exponential suppression 
of the amplitude of long orbits in the semiclassical Green function. The special half integer decays can 
probably be always associated with the presence of impurity and wedge diffraction. 

8 Summary 

In this paper we analyzed the conductance of mesoscopic devices where diffraction plays a major role. 
Such pseudo integrable systems are typical, when the form of the device is human designed. We derived 
formulas expressing the reflection and transition matrix elements for one and many diffractive points. 
We verified this for the special case of impurity and wedge diffraction. Diffraction can cause back 
scattering in situations, where usual semiclassical back scattering is absent causing an erodation of ideal 
conductance steps. The length of diffractive periodic orbits and diffractive loops can be detected in the 
power spectrum of the reflection matrix elements. The large I part of the power spectrum shows l// 1 / 2 
decay due to impurity scattering and 1 / Z 3 / 2 decay due to wedge scattering which we think is a universal 
sign of the presence of diffractive scattering in pseudo integrable waveguides. 
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The probability to hit a point like impurity is zero. 

The regularization removes the zero length trajectory from the point to itself, which causes the 
divergence of the Green function otherwise. 



